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ASTROPHYSICAL AND GEOPHYSICAL VARIATIONS IN C 14 PRODUCTION 
I. INTRODUCTION 

We shall study time dependent variations in the pro- 
duction rate of C ^ resulting from changes in various astro- 
physical and geophysical factors* The general success of the 
C 14 dating method implies that to a first approximation the 
production rate of C ' 4 has been essentially a constant for the 4 
last several millenia. However, a variety of phenomena cause 
changes in the production rate which in turn can produce 
measurable perturbations in the biospheric C ^ 4 activity. In 
particular we shall consider changes in production caused 
by var i at i ons ' i n the terrestrial and interplanetary magnetic 
fields leading to changes in the flux of galactic cosmic rays 
at the earth; by enhanced fluxes of particles from solar flares; 
and by variations in the local interstellar cosmic ray flux 
produced by nearby supernova explosions. 

We calculate the yields of C ^ 4 from the interaction of 
protons, alpha-particles and high-energy gamma rays as functions 
of energy. Using models of the modulation of the cosmic ray 
spectrum by terrestrial and interplanetary magnetic fields and 
the possible time variations of these fields, we then compute 
the resultant changes in C * 4 production. Similarly, from 
measurements of the solar flare particle fluxes, and a diffusion 
model for possible cosmic ray variations from nearby supernova 



sources, we calculate the possible production variations. 

Finally, we calculate the resultant changes in the biospheric 

1 k 

C activity using a simple, two- reservoi r model for the dis- 
tribution of exchangeable carbon. From a comparison of these 
calculations wi th the observed fluctuations of activity, 
we attempt to deduce the most plausible causes of these 
fluctuations and at the same time place upper limits on the 
variations of the terrestrial and interplanetary fields and 

the intensities of impulsive cosmic ray sources within the 
* , c 

last 10 years. 

II. C 14 PRODUCTION 

In order to study variations we must first determine 
the dependence of the production rate on the cosmic- ray 

, 1 L 

energy spectrum. The rate of C production by cosmic rays 
incident on the atmosphere may be written: 

__ 00 

Q=Vj TTcp.(E)m.(E)dE , (l 

i ° 

where Q is the production rate in atoms per second in 
a cm column of atmosphere; i is the species of cosmic- 
ray particle (in this calculation only protons, He, C, N, 
and 0 nuclei are important); cp - ( E ) is the differential flux 
of the i*th particles per cm • sec ster Mev per nucleon Inci- 
dent on the atmosphere; and mj(E) is the total production or 
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yield of resulting from the interaction in the atisios- 
phere of a particle of species, I, and initial energy per 
nucleon, E. 

Considering for the sake of brevity only the first and 

* 

second generation interactions In the bMild up of the nujcjepn—. 

~ *. • « • * * 

• • 

- ~ ** ^•■"CTTscSTO’ thfs~yield may be written ( L i ngenf el ter et al . 1965): 
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The first integral term in the equation is the yield of 
produced by the capture of first generation neutrons 
produced by the incident particle, i. The second term is 
the yield of produced by the capture of second genera* 
tion neutrons, produced by all first generation nucleons, 
s, which were produced by the incident particle, i. Similar 
more lengthy terms describe the production from subsequent 
generations. P - ( E ) is the probability that a neutron, 
resulting from interactions of the incident particle i of 
energy E and its secondaries, will be captured to produce 
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C 14 . Diffusion calculations by Lingenfelter and Flamm (1964) 
show that this is a rather slowly varying function of the 
incident proton energy E varying from 0.35 at 100 Mev and 
below to 0.65 at a few Bev and above, cr. and cr are the 

I D 

cross sections for neutron production by the incident particle, 
i, and the secondary nucleons, s, and a. is the cross 

I D 

section for secondary nucleon production by the incident 
particle, f. (E ,E') is the probability that the secondary 

I 3 3 

nucleon produced by the incident particle of energy E' will 
have an energy between E and E + dE . R(E) and dE/dr are 

ta) 3 J 

the range and differential energy loss of the particles in 

> 

the atmosphere. L( E ' , E ) is the mean attenuation length of 
the particles in the energy interval E' to E. 

Data on the cross sections for neutron production by 
protons in air are presented by Lingenfelter et al . (1965). 

For lack of data on the neutron production in interactions of 
cosmic- ray He and heavier nuclei with air we have assumed 
these cross sections to be A a times the cross section for 
proton- i nduced production at the same energy per nucleon, 
where A is the atomic number of the cosmic ray nucleus and 
a ranges from 2/3 at energies less than 100 Mev/nucleon to 
1 at energies greater than 1 Bev/nucleon. We have further 
assumed that the total cross section for secondary nucleon 
production is twice that for neutron pv eduction. The data 
on the other parameters have been summarized by Lingenfelter 
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and Flamm (1964) and Lingenfelter et al * (1965) 

The resultant yields per particle, calculated from 
equation (2), are shown in Figure 1 as a function of the inci- 
dent particle rigidity for cosmic-ray protons, He and CNO 
nuclei incident on the atmosphere. The first generation 
neutrons account for most of the production for incident 
proton energies less than 4 Bv. Subsequent generations con- 
tribute 0.1 of the total production at 1.5 Bv and 0.5 at 
4 Bv. At rigidities greater than 10 Bv the yield is 
assumed to be proportional to the rigidity to the 0.4 power. 

Using the measured energy spectra and nuclear abundances 
of the cosmic ray flux at the top of the atmosphere at solar 
minimum and solar maximum (Webber, 1967) we can evaluate 
equation (1) for the polar production rate at both times. 

The calculated solar maximum to minimum production ratio 
of 0.67 is in excellent agreement with the value of 0.70 
determined by Lingenfelter (1963) from atmospheric neutron 
measurements at solar minimum and maximum. Furthermore, 
using vertical cut-off rigidities as a function of geomagnetic 
latitude based on a dipole filed (p^(Gv)= 14.9 cos\), we 
can calculate the latitude dependence of the production 

Q X =I J E . TTcpj(e)mj(E)dE . (3) 

i ix 

Where is the cut off energy corresponding to p^ . 
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Because of the different rigidity dependence of the yield 
for protons and multiply changed nuclei (see Figure 1) the 
contribution of the latter to the total production varies 
from 32% of the poles to ^8% at the equator. The relative 
latitude variations of the total production thus calculated 
are in very good agreement with those determined by Lingenfel ter 
(1963). However, because of the compounding of the uncertainties 
in the cross sections and other parameters, when the inter- 
actions of second and higher generations of nucleons become 
important at energies greater than a few Bev, it is necessary 
to normalize the absolute yields at high energies in order 
to give agreement .wi th the best estimate of the production 
rate. 

The previous determination of the production rate 
by Lingenfel ter (1963) was based on normalization of the 
relative production rate to atmospheric thermal neutron 
measurements. In particular, the altitude and latitude 
dependence of the cosmic- ray neutron production was determined 
from data on the rates of cosmic- ray nuclear-star production 
and ionization, and on the equilibrium flux of cosmic-ray 
neutrons in the atmosphere. From this source distribution, 
arbitrarily normalized to a production rate of 1 neutron per 
cm column of atmosphere per second at the geomagnetic pole 
during solar minimum, the energy-dependent, equilibrium 
neutron flux was computed using the diffusion equation. The 
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relative altitude and latitude dependence of the pro- 
duction rate was then calculated from this flux distribution. 

The absolute normalization of the production rate was de- 
termined by normalizing the calculated thermal neutron flux 
to experimentally determined rates, measured by balloon-born, 
bare counters. Since most of the measurements were 

not made at solar minimum on maximum, it was arbitrarily 
assumed that the variation of the relative flux between times 
of solar maximum and minimum was linearly proportional to 
the sunspot number. Now, however, it is possible to directly 
determine the relative flux as a function of time from the 
equations given above and the measured variation of the cosmic- 
ray intensity over the last solar cycle. 

Lockwood and Webber (1967) have shown that the polar cosmic- 
ray flux variation between solar minimum and maximum can be 
wri tten 

cp(P,t) - cp(P,t 0 )exp{n(t)/Pa} , (4) 

where P and p are the cosmic- ray particle rigidity and 
relative velocity, respectively; tp is the time of solar 
minimum; and r\(t) is a function ranging from 0 at tp to 
about -2Bv at solar maximum. The solar modulation function, 
n, can be directly related to the corrected daily mean Mt. 
Washington neutron monitor counting rate (Lockwood, personal 
communication) and Cheltenham ionization chamber rate 
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(Forbusb, 1958); the combined record of these monitors covers 
a period of three solar cycles back to 1937. Using this tempo- 
ral dependence of the flux, we have calculated the relative 
variation of the equilibrium thermal neutron flux which is 
directly proportioned to the production rate defined in 
equations (2) and (3). From this variation we can then 
renormalize the relative neutron and production rates of 
Lingenfelter ( 1963) . 

For this renormal izat i on we shall use the thermal 
neutron measurements in the region of the Pfotzer maximum 
by Yuan (1951) on June 9, July 25, 1948, and January 8, 1949 
at, 51 .8° N. Geomag. Lat.; Soberman (1956) on August 27, 29 
1952 at 88.6° N. Geomag. Lat., August 24, 28 1954 at 55° 

N. Geornag Lat. and September 6, 8, 9, 1953 at 10.1°; Reidy 
et al . (1962) on August 23, I960 at 49° N. Geomag. Lat., 

Smith et al . (1962) on July 20, 1961 at 57° N. Geomag. Lat»; 
and Miles (1964) on February 6, July 5, 6 and October 22, 23 
1962 at 41° N. Geomag. Lat. We have not included the air- 
craft-born measurements by Hess et al . (1959) which appear 

to have been affected by local neutron production in the air- 

2 

craft. The absolute neutron production per second in a cm 
column of atmosphere at the geomagnetic pole during solar 
minimum determined from the normalization of each set of 
measurements is: 1.1 ± 1.6 for Yuan (1951); 9.4 ± 1.4 for 
Soberman (1956); 8.3 ± 2.1 for Reidy et al . (1962), 9.7 ± 2 
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for Smith et al . (1962) and 7.2 ± 1.4 for Miles O a 64). The 
weighted mean of these values is 8.35 ± 1.7. This is 7% 
lower than the value of 9.0 ± 1.8 previously determined by 
L i ngenf el ter ( 1963 ) . 

This renormalization gives a global average production 

rate of 2.42 ± 0.48 C^/crn^ sec during solar minimum (1953-4) 

and 1.93 ± 0.39 during solar maximum (1957-8). The latitude 

dependence of the neutron and production rates dur ing 

solar minimum and maximum are tabulated in Table 1, The 

dependence of the global average production rate on the 

solar modulation function, r\ , and hence on the daily mean 

neutron monitor and ionization chamber rates is shown in 

Figure 2. Using quarterly average values of the Mt. Washington 

daily mean neutron monitor rate we calculate that the average 
1 4 

C production over the last solar cycle (1955 to 1964) was 
• 1 4 2 

2.1 ± 0.4 C /cm sec. Including the Cheltenham ionization 
chamber record back to 1937, we find a three cycle average 
(1937 - 1967) of 2.2 ± 0.4 C ,4 /cm 2 sec. 

III. C 14 ACTIVITY 

In order to investigate variations in the c' 4 activity, 
it is necessary to assume a model for the distribution of 
exchangeable carbon in the terrestrial reservoir. The 
relatively long term variations which we discuss quantita- 
tively in the present paper can be investigated by using a 
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two reservoir model: an upper reservoir containing A grams 
of carbon per cm of the earth's surface, which represents 
the atmosphere and biosphere and possibly the upper hundred 
meters of ocean above the thermocline; and a lower reservoir, 
containing B grams of carbon per cm , which represents the 
oceans (Wood and Libby, 1964; Ramaty, 1965; Houtermans, 1966). 
The ratio B/A is defined as a parameter v . If the upper 
hundred meters of ocean are part of the upper reservoir v«30, 
whereas if the entire oceans belong to the lower reservoir 
v«6G . 

Let R a and R g be the radiocarbon concentrations in the 
upper and lower reservoirs, respectively. The differential 
equations satisfied by the contents AR^ and BR^ are (e.g. 
Ramaty, 1965) 


m < ar a> = 1 • : " A 


AR, 

T 


AR, 


BR 
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(5) 
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(br b ) = - 
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where: Q is the global average C production which may be 
time dependent; t is the mean life of radiocarbon equal to 
8300 years; T is the average residence time of a carbon atom 
in the upper reservoir equal to about 25 years (Wood and 
Libby, 1964); and S is a possible loss rate from the lower 
reservoir due to sedimentation. 


( 6 ) 
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By assuming that S is time independent and by defining 
the decay rates in the upper, lower and total reservoirs, 
respectively, as 

J A = \ (A+B)R a + S(1-T*/t) 


J B = \ (A+B)R b + S( )+T*/vt) 


J t (AR a + BR b ) + s 


where 


T * - / 1 . 1 . 1 X 

T - + t + vr 


-1 


the solutions of equations (5) and (6) can be written as 

J A (t) = i j" e (t '- t)/T Q(t')dt' + | f e (t '- t)/T *Q(t')dt' 


J B (t) = -l J " t ^ T Q(t')dt' - 1 J " t ) /T Q( t ‘ )dt 

-00 -00 

> 

J t (t) =1 J l e (t, ' t)/T Q(t')dt' 

— 00 

In order to evaluate J^, J g or J t 'at a given time t, 
it is necessary to perform the integrals in equations (10), 
(11) and (12), and, therefore, Q(t') must be known over a 
time period of at least several mean lives t prior to t. 

Since the geophysical and astrophysi cal factors which 


(7) 

( 8 ) 

(9) 

( 10 ) 

(ID 

( 12 ) 
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determine the C production are known at best only at the 
present and possibly over a relatively short time period in 
the past, arbitrary assumptions would have to be made about 
the radiocarbon production in the more distant past. How- 
ever, we can eliminate the need for such arbitrary assumptions 
by introducting a parameter, R Q , equal to the present ratio 
between production and decay in the total reservoir, 


R 0 « Q(0)/J t (0) , (13) 

and by transforming the time variable, t, into a new variable, 

u, which has its origin at the present and increases backward 

in time. If Q is constant for several mean lives, t, an 

equilibrium is established between production and decay, 

1 4 

and R q = 1 . The general success of the C dating method 
would imply that even if Q is somewhat variable, R Q must be 
close to unity. 

In terms of the new time variable, u , and the parameter, 

R q equations (10) and (11) can be solved for J^(u)/J^(0) and 
J b (u)/J b (0): 


J A,B (u) _ , 11 ) _ 1 r „ , u-u\ Q( 

p( ") -{ exp( — 


U 


q( ) 
H5T 


du' - 


B 


-7 <'•"'>[/ 

u+U 

-I e * P <T^ du '] 


Q( du' 




(14) 
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where the ratios Q(u)/J^ B (C) are given by 


Q(u) _ Q(u) 

~wn 


1 



exp( - 



Q(u') 

Tjnrr 


du' 


( 15 ) 


and the parameters v and -1 correspond to A and B, respectively. 

U tends to infinity, but to make the result exact it is 

•ft 

sufficient that it be only much larger than T . Thus for a 

given R q , equations (14) and (15) can be evaluated at any 

time u in the past provided that the production function Q(u') 

is known over the time interval starting at the present and 

extending somewhat beyond the time considered. The unknown 

C 14 production at earlier times is represented by R q . 

The allowed range of values for the parameter, R , can 

be determined from direct C^ 4 measurements and the production 

calculations discussed above. According to Karlen et al . (1966), 

the specific activity of 19th century wood' (which in our model 

equals R^/t) is 13-56 ± 0.07 disintegrations per minute per 

gram of carbon. The sedimentation rate S was estimated by 

1 4 2 

Libby (’965) to be 0.5 ± 0.3 G atoms per cm second. Assuming 
that A+B = 8.3 g/cm (Libby, 1965), from equation ( 7 ) we get 
- 2.38 ± 0.3 disintegrations per cm second. Using the 
average Q - 2.2 ± 0.4 determined above, we find that the present 
value of the ratio Q/J^ (normalized to 19th century conditions) 
equals 0.96 ± 0.29- Solving for R Q , from equation (15) we obtain 
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Q(Q)/J AiB (Q) 




(V.-I) 



where 
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Since Q(u')/Q(0) is not expected to vary very much over 
several mixing times T , we take T = T =25 years. For 
v=60 , t = 8300 years, and Q(0)/J^(0) = 0.96 ± 0.29, we 
find that 0.75 £ R Q £ 1 .61 . This range of values of R Q 
reflects the uncertainties in the C^ 4 production, Q, and 
removal by sedimentation, S. Since these uncertainties are 
large, they introduce in turn large uncertainties in the 
allowed values of R Q . However, as mentioned above, the 
relative constancy of the C^ 4 activity in the past suggests 
that R q must be confined to a much narrower range of values. 
We shall assume, therefore, that R Q is essentially a free 
parameter to be determined from the study of the long term 
C 14 activity variations. 


IV. c' ^ VARIATIONS 

We shall consider now, the various geophysical and astro- 

1 h 

physical factors that may influence the global average C 
production. Since radiocarbon production and decay appear to 
be close to equilibrium, we can first treat those factors 


( 16 ) 


(17) 
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which may modulate a constant cosmic ray background and then 
consider separately the variations in such a background aris- 
ing from discrete cosmic ray sources. The modulations which 
we shall consider are those associated with variations in the 
geomagnetic and interplanetary magnetic fields, and the dis- 
crete sources which we shall treat are solar flares and super 
nova explosions. 

Geomagnetic and Solar Variations 

The effect of a varying geomagnetic field on pro- 
duction was considered by a number of authors (E! sasser et 
al . 1956; Ramaty, 1965 Wada and Inous, 1966,* Kogoshi and 
Hasegawa, 1966; Bucha and Neustupn^, 1967). Basically the 
production at a given geomagnetic latitude \ depends on the 

magnitude of the vertical cutoff rigidity at that latitude, 

2 k 

P(\) = M/4r cos 1 , where r is the earth’s radius and M is 
its dipole moment. Because of the rapid mixing of radio- 
carbon over the earth’s surface and because of the isotropy 

of the cosmic rays, the only variations of interest are in 

* 

the magnitude of the geomagnetic dipole and riot of its 
d i rect ion . 

Variations in the geomagnetic dipole moment over the 
last 9000 years have been recently summarized by Cox (1969). 
The reported values of M represent averages of measurements 
made at different locations under the assumption that the 
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geomagnetic field is a perfect dipole. Since this assumption 
is known to be only approximately valid, it will introduce 
some uncertainty in the deduced variations. Furthermore, 
many of the samples used for the geomagnetic measurements 
were dated by radiocarbon, which, because of the geomagnetic 
variations themselves, may introduce additional uncertainties 
in the deduced variations, nevertheless, lacking more 
accurate geomagnetic data, we shall use the data as summarized 
by Cox (1969) (see Table 2). 

In addition to geomagnetic variations, the cosmic ray 
flux incident on the top of the atmosphere is also influenced 
by the state of the interplanetary magnetic field. As shown 
in equation (4) above, this can be characterized by a para- 
meter n , which had the values 0 and -2GV for solar minimum 
and solar maximum, respectively, of the last solar cycle. From 
Figure 2 the mean production over three solar cycles of 
2.2 atoms pe r cm second corresponds to r\ = - .7GV. In the 
absence of direct measurements outside the solar system, the 
value of r\ that characterizes the cosmic ray flux in inter- 
stellar space can only be estimated by using indirect arguments. 
By comparing the measured fluxes of deuterons and helium-3 
nuclei in the cosmic rays with those expected from nuclear 
interactions in interstellar space, Ramaty and L i ngenf el ter 
(1969) found that the mean cosmic ray flux in the galaxy can 
be related to that observed at the earth by a modulating 
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parameter r\ * 0.35 ± 0,15 GV. This value is consistent with 
the cosmic ray gradient measurements of 0‘Gallagher and 
Simpson (1967), who found a An of 0.2 GV between the orbits 
of the Earth and Mars. 

The C 14 production rates as functions of M, for r\ equal 
to -2GV, -.7GV, 0, .35GV and ,5GV, are shown in Figure 3» 

For values of M close to its present value (8x10^ gauss cm^) 
the production varies approximately like M’ 0 *^ which is in 
agreement with results previously obtained by El sasser et al . 
(1956) and Ramaty (1965). However, as is evident in Figure 3, 
this simple relationship breaks down for much smaller and much 
1 arger val ues of M. 

Suess (personal communication) has suggested that, since 
the cutoff rigidities at all latitudes decrease as M decreases, 
the global average C^ 4 production becomes more sensitive to 
solar modulation variat ions. From Figure 3 we see that the 
change in the global average production rates from solar 
minimum (0GV) to maximum (-2GV) would have varied from 
AQ = 0.48 C^ 4 cm ^ sec ^ at about 1500 B.P. , when the dipole 
moment had a maximum value of about 11. 4x1 gauss cm^ , to 
AQ = 0.71 C^ 4 cm ^ sec”^ at about 5500 B.P., corresponding to 
a minimum value of 5.1x10^ gauss cm^ . Although the eleven 
year solar modulation effect on C^ 4 activity may be obscured 
by solar-flare produced increases in C 14 production, as we 
shall discuss later, both the solar flares and longer period 
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solar modulation variations would produce changes in C^ 4 
activity, the magnitude of which would also depend on dipole 
field variations. 

Since there is no information on the long-term variation 

1 f 

of r) , we shall consider only the C activity variations 
that may result from the measured geomagnetic variations. 

Using the dipole moment variations (Cox 1969), listed in 
Table 2, and the resultant production variations for the mean 
value of rj ( - . 7GV) , shown in Figure 3, we have evaluated 
equations (14) and (15) for v=60, T =25 years, and a range 
of values of R q . As mentioned above, for a given value of 
R q the C 14 activity variations can be uniquely determined over 
almost the entire time period for which the magnetic data is 
available. The resultant C^ 4 variations in the upper 
reservoir are shown in Figure 4 for R Q = 1 and 1.05 together 
with the measured radiocarbon variations based on dendrochrono 
logical studies (Suess, 1965, 1967). 

The minima and maxima at about 1000B.P. and 5500 B.P. 
correspond, respectively, to the maximal and minimal dipole 
moments at about 1500 B.P. and 6000 B.P. The magnitude of 
these extrema depend critically on the value of R Q . As can 
be seen from Figure 4, the measurements are reasonably well 
bracketed by 1 .0 £ R Q < 1.05 and thus the gross features of 
the C 14 activity variations can be understood in terms of the 
variations of the dipole moment. 
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The limits thus placed on R Q contain information on 

both the history of cosmic ray intensity and on the rate of 

radiocarbon sedimentation. For this variation of the dipole 

moment, a cosmic ray intensity outside the geomagnetic field 

that is constant in time would give an R Q equal h to ".04. 

Lower values of R q imply a higher cosmic ray intensify 

time in the past and higher values imply a lower intensify. 

We shall consider this point in the later discussion of possible 
14 

C variations resulting from supernova explosions, and turn 
now to the implications for the sedimentation rate. 

For the upper limit R Q < 1.05, we find from equation (16) 

that 

Q(0)/J a ( 0) = R 0 /(l + ^L- R q ) < 0.88 

and since, as determined above, Q(0) ^ 1.8 C 1If cm 2 sec ‘, 

J A (0) £ 2.04 C 1 ** cm” 2 sec"*. This is consistent with the 
range 2.38 ± 0.3 C 1 ** cm” 2 sec/”* given above. However, if 
removal by sedimentation is neglected, J A 1.88 C 1 ** cm 2 sec * 
and this is obviously Inconsistent with the lower limit on J A 
of 2.04 C 1 * cm” 2 sec” * obtained from the upper limit on R Q . 
These arguments would require, therefore, that radioactive 
carbon be removed from the exchangeable reservoir at a rate 
of at least 0.16 C llf cm' 2 sec" 1 . From the argument presented 
above, however, it would appear that by decreasing the product 
vt* one would increase Q /J A and therefore decrease the lower 


-» Tv 

limit on , A much lower value of vT however, would be 

1 L 

inconsistent with the C measurements in the deep oceans. 

In fact, by solving for J A (0) and Jg(0) in equation (16), we 


obta ? n 


J B<°> _ '-f- R o 




't A, 


which, evaluated for v=60, T = 25 years, and R q = 1.025 ± 0.025, 
gives Jg/J^ ® 0.84 ± 0.005. This is consistent with the measure- 
ments of Bien et al . (1962) who found that the radiocarbon 
activity in the deep oceans is lower than that of the biosphere 

' — it 

by about 15% to 23%. If, however, vT were lowered by a factor 
of about 2 (which would be required for consistency without 
sedimentation), Jg/J^ would equal 0.92 and this is in dis- 
agreement with the C 14 measurements in the deep oceans. 


Solar Flare Variations 

The production of C ^ 4 by the interaction of solar-flare 
particles with the earth's atmosphere has been studied by 
Simpson (I960), Lai and Peters (1962) and Lingenfelter and 
Flamm (1964). The spectrum of solar-flare particles arriv- 
ing at the earth has been measured by Freier and Webber (1963) 
who have shown that it can be represented by an exponential 
in particle rigidity, P. 


cp(P) = cp Q exp{ -P/P q} 
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The characteristic rigidity, Pq , varies with different flares. 
For the major events of the last solar cycle it has ranged 
from 50 Mv to 325 Mv. Measurements of the solar particle 
abundances by Freier (1963) and Biswas et al . (1963) have 
also shown that the ratios of the differential fluxes per 
unit rigidity of protons, alpha particle and medium nuclei 
are essentially independent of rigidity. However, since the 
solar particles are predominatly of rigidities less than 1 Bv, 
it can be seen from Figure 1 that C^ 4 production by He and 
CNO is negligible compared to that by protons. 

Using the above form of the particle spectra and integrat- 
ing equation (3) over dcosX , we calculate the global average 
C^ 4 production rate. Moreover, during geomagnetic storms, 
when the solar flare particles arrive, the vertical cut off 
rigidities as a function of latitude are different from those 
of the average undistributed field seen by the cosmic rays. 
Geomagnetic storms following major solar flares can reduce 
the effective cutoff rigidities to as small as 20% of the 
normal cut off rigidities. Thus to place bounds on the solar- 
flare C 14 production we have calculated the production rate 
for both the normal dipole cut off rigidities and for reduced 

cut off rigidities. These production rates, normalized to an 

- 2 - 1 

incident solar particle flux of Ip cm sec of energy greater 

than 30 Mev, are shown in Figure 5 as a function of the charac- 

1 4 

teristic rigidity, Pq . As can be seen, the C production 
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rate is much more st rongl y .dependent on the characteristic 

rigidity than on the cut off rigidity. 

Malitson and Webber (1963) have surveyed the particle 

intensity measurements for the major solar flares of the 

last solar cycle and have estimated both the integral flux 

at the earth and characteristic rigidity for each event. 

Using these values, listed in Table 3, we have calculated 

1 4 

for each flare the gl obal -average C production, averaged 

over a year for comparison with the cosmic ray production 

1 4 

rate. The total C production per 1 unit surface area for any 
event is thus 3x1 0^ times the value tabulated. The solar- 
flare C^ 4 production rate averaged over the 11-year solar 
cycle lies between 0.31 and 0.12 C^ 4 crrf^ sec“* depending 
on the reduction of the cut off rigidity. This is between 
14% and 6% of the solar-cycle averaged cosmic- ray C^ 4 pro- 
duction and as can be seen the bulk of the C^ 4 is made in 
one or two events. Moreover such solar flares produce enough 
C 14 to greatly modify the solar-cycle dependence of the total 
C 14 production rate (see Figure 6). Therefore, for particularly 
active solar cycles, such as the last, the 11-year periodicity 
resulting from solar modulation is compl etel y obscured by solar 
flare effects. As a result for very active cycles one should 
not expect to find anti-correct ion between the atmospheric 
C 14 production and sunspot numbers or other indicators of the 
11-year solar-cycle activity variation. However, this does 
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not preclude correlations with longe** period variations in 
solar activity, such as those suggested by Stuiver (1961). 

The increase in atmospheric C^ 4 activity AJ(t) resulting 
from an impulsive incremental increase in C 14 production AQ 
at a time t=0 may be written 

AJ A (t) AQ v _„_ r WT *, 

_j_ « tj- - exp{-t/T } 

where T « t and and Q are time averaged values. Using 
the values of these quantities, defined as above, we see 
that the solar flare of 23 Feb. 1956 may have produced as 
much as 0.75% increase in the specific activity of atmospheric 
C 14 and all flares of the last solar cycle may have produced 
as much as a 1.1% increase. 

Thus measurable solar-flare produced increases in 
1 4 

atmospheric C activity may be associated with periods of 
great solar activity in the past and such a source may account 
for some of shorter time variations observed in C^ 4 activity 
in dendrochronol ogi cal 1 y dated samples. Unfortunately an 
experimental measurement of solar-flare produced increases 
during the last solar cycle is not possible because of the 
much larger increases in atmospheric C^ 4 activity produced 
by nuclear testing during the same period. However, such 
measurements are possible for known, earlier flares, although 
the measurements of the energetic particle intensities and 


( 20 ) 
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spectra ore not available for them, 

Supernova Vo r i at i ons 

Lastly we shall consider possible C^ f variations caused 
by supernovae which collectively are the most likely source 
of most cosmic ray particles. There are two types of increases 

i / 4 

in C‘ production which may result from a relatively nearby 
supernova explosion: a short-term increase produced by a 
possible gamma ray burst associated with the explosion, and 
a much longer term increase and subsequent decrease resulting 
from enhancement of the local background cosmic ray flux by 
the arrival of cosmic rays accelerated in the explosion. 

The possible C^ 4 increase produced by a supernova gamma 

ray pulse has been discussed in detail by Konstantinov and 

Kocharov (1965, 1967). They point out that gamma rays will 

interact in the atmosphere and through photonuclear reactions 

produce neutrons which in turn will be captured by nitrogen 

to yield C 14 . Although the photoneutron cross sections for 

nitrogen and oxygen show a peak at about 25 Mev and decrease 

at higher energy, the development of a photoel ect ron-brems- 

strahlung cascade causes the neutron and hence C [ yield to 

3 1 4 

be essentially constant at about 10 C per erg of gamma 
rays (> 10 Mev) incident at the top of the atmosphere. This 
yield in terms of C atoms per gamma ray photon is comparable 
to that for cosmic ray protons of the same energy. Estimates 
of the total energy emitted in gamma rays during a s up r nova 
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explosion range from 10 4 ® ergs (Colgate, 1968) to 10^ ergs 
(Gould and Burbidge 1965). Konstantinov and Kocharov (1967) 
have shown that measurable increases in the atmospheric C^ 4 
activity would have been produced by historical supernovae if 

hq 

their gamma ray emission energies were greater than 10 y ergs 
(see Table 4). Like the C^ 4 increases associated with solar 
flares, these increases would decay with a mean life of T 
At present no systematic search has been made for such 
increases, but such a search would prove quite valuable in 
at least setting an upper limit on the supernova energy 
emitted in gamma rays. 

We now turn to perhaps the most fundamental variation in 
C 14 activity, that reflecting local variations in the cosmic 
ray intensity resulting from nearby supernova sources. In 
general the cosmic rays may be assumed to come from many super- 
nova randomly distributed throughout the galactic disk. Studies 
of the production of secondary isotopes, such as D and He - ^ 
(Ramaty and L i ngenfel ter , 1969) and the light elements (Shapiro 
and Silverberg, 1968), by cosmic ray interaction in the 
interstellar medium have shown that the mean amount of matter 

_ o 

traversed by local cosmic rays has been between 3 and 4 gm cm 

In the interstellar gas of the galactic disk, which has a 

hydrogen density of the order of 1 atom cm , this path length 

6 

corresponds to a mean lift of about 2x10 years for relativistic 
cosmic rays. If we assume that this is the mean life time t 
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for cosmic ray escape from the galactic disk, then the 
equilibrium cosmic ray energy density, w, may simply be 
wri tten 


w 


fW SN T e 


» 


where f is the frequency of supernova explosions in the 

disk; is the total energy of cosmic rays accelerated 

in a supernova explosion and V is the volume of the disk 

(4x10 crrr). The frequency of supernova explosions in the 

galaxy has been estimated to be about 2x10 per year 

(Katgert and Oort, 1967, Kesteven 1968). Therefore the local 

- 1 2 - 3 

cosmic ray energy density of 10 ergs cm requires that 

the total energy of cosmic rays accelerated in a supernova 
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explosion, be of the order of 10 ergs. 

Assuming isotropic, three dimensional diffusion, the 
cosmic ray flux cp(E,t,r) arriving from a source of age, t, at 
a distance, r, may be written: 

<p(E,t,r) = W $N ^LEI exp{-r 2 4 Xvt} . 


where N(E) is the number distribution of cosmic rays at the 
source per erg of total cosmic ray energy; v is the cosmic 
ray particle velocity; and X is the cosmic ray diffusion mean 
free path in interstellar space. This expression may be re- 
written in terms of a dimensionless time, t/t , where t is 

m m 
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the time to maximum for a relativistic part icl e, equal to 

r /2\c . Assuming that N(E) ~ (me + E) , we get cp 

2 

(particles per m sec»ster*Hev per nucleon), 

W CKI ( 10^ erg) o c/h t -,3/2 

cp-l6-§S — p(1-p 2 ) 5/4 [-2 e X pt- -2}] . ( 

r^( 1 00 pc) L Pt 0t J 

Substituting this flux dependence into equations (2) and (3), 

we then calculate the global average production resulting 

from such a source. This shown in Figure 7, normalized to 
SO 

W$n = 10 ergs and r=l 00 pc. 

The supernova observed in historical times (listed in 
Table 4) are all so distant that their age is much less than 
the diffusion time (t/t « 1) and therefore cosmic rays from 
them have not yet reached the earth. However, the recent 
discovery of pulsars, which may be neutron-star remnants of 
supernovae (Gold, 1968 and Pacini, 1968), possibly provides 
new clues to the age and distance of nearby supernova remnants. 
Grewing and Priester (1969) have estimated the age of pulsars 
from the observed periods and the rates of slowing down and the 
distance can be estimated from the dispersion measure. Survey- 
ing the available pulsar data (Maran and Cameron, 1969), we 
find that for a constant the largest cosmic ray flux at 
present would be coming for PSR1929+10. This pulsar has a 
period of 0.227 sec, which from Grewing and Priester's studies 
would suggest an age of about 10^ years, and a dispersion 
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measure of 8 electrons pocrrf^, which would give a distance 
of 80 pc for a free electron density of 0,1 cm~^. These 
values give a t/t of 10 at the present. The maximum flux 
at the earth from this source would have arrived about 9xl0 4 
years ago and the flux would now be decreasing with the -3/2 
power of the age. 

Over the last 8000 years the flux from PSR1929+10 would 

have decreased by 13%. If we assume that this source is the 

principal perturbation on more constant background arising 

for all more distant sources, then the decrease in the total 

cosmic ray flux during this period would depend on the fraction 

of the total flux presently contributed by this single source; 

this in turn depends on the total cosmic energy emitted by 

the source. As discussed above the values of R which are 

o 

consistent with the dend rochonol og i cal 1 y determined C 14 varia- 
tions, suggest that the C^ 4 production rate 8000 years ago 
could have been as much as 5% greater. A nearby source such 
as P S R 1 9 29 + 10 could produce increases of this magnitude. For 
example W<^= 3.76x10^ ergs could produce an increase of 5% 

C A 

and a lesser increase of 1% would result the W SN = 0.75x1 0 3 ergs. 
This can be seen in Fig. 8, which gives the time history of the 
C 14 production rate for a local cosmic ray source of the 
distance and age of PSR1929+10 superimposed upon a constant 
background. If the cosmic ray flux was as much as 5% higher 
8000 years ago and if this excess resulted solely from a 10^ 


year old nearby source, then the cosmic rays from this source 
would make up 38% of the total flux at the present. If the 
increase were only 1% then the source would only contribute 
only 8% of the present flux. It is of course much more likely 
that, if such an increase does exist, it is due to more than 
one source and from the same diffusion model we see that 
increases of 5% to 1% would also be expected if all of the 
present cosmic ray flux came from sources that were 0.25 to 
1x10 years old, 

V. SUMMARY 

In summary then we have renormalized the absolute C^ 4 
production rate based on new atmospheric neutron measure- 
ments and from the cosmic ray neutron and ionization monitor 
records, we have determined the average production rate to 
be 2.2 ± 0.4 C 14 cm ^sec ^ for the last three solar cycles 
(1937 to 1967). Because of the uncertainties in the calcula- 
tion of both the production rate and decay rate of C 14 we 
find that the best determination of the ratio of these two 
rates is obtained from the C^ 4 variations determined from 
dendrochronology. We have shown that the major component of 
these variations can be understood in terms of measured geo- 
magnetic field variations over the last 1 0 4 years. Shorter 
time variations may result from solar modulated cosmic rays, 
solar flares particles and possibly supernova gamma rays, while 
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longer time variations may reflect changes in the local 
cosmic ray flux resulting from nearby supernova explosions. 
All of these variations have been treated quantitatively and 
further measurements of variations can give valuable 
information on these processes. 
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TABLE 1 . 


Cosmic-Ray Neutron and C 
(n and C ^ cm"^ sec - ^ ) 


14 


Production 


Solar Minimum, 1953 

-1954 

Solar Maximum, 1957-1958 

Geomagnet i c 
Lati tude 

Neutron 

Source 

Carbon 14 
Product ion 

Neutron 

Source 

Carbon 14 
Production 

0° 

1 .37 

0.91 

1 .31 

0.86 

10° 

1 .42 

0.94 

1 .35 

0.89 

20° 

1 .72 

1.13 

1 .6) 

1 .07 

30° 

2.56 

1 .70 

2.28 

1 .51 

40° 

4.27 

2.80 

3 .46 

2.27 

50° 

6.53 

4.20 

4.95 

3.19 

60 ° 

7.92 

4.88 

5.61 

3.50 

70° -90° 

8.35 

4.99 

5.61 

3.50 

Global average 3.81 

2.42 

2.99 

1 .93 





TABLE 2. Geomagnetic Field Variation 

(Cox, 1969) 


Years 
B. P. 

Dipole Moment 

25 3 

10 gauss crrr 

Years 
B. P. 

Dipole Moment 

25 5 

10 gauss cm 

0 

8,0 

5500 

5.1 

500 

10.5 

6000 

5.1 

1000 

11.2 

6500 

6.4 

1500 

1 1 .4 

7000 

8,0 

2000 

11.3 

7500 

9.5 

2500 

10.2 

8000 

10.2 

3000 

9.3 

8500 

1 1 .2 

3500 

9.0 

9000 

11.5 

4000 

8.7 

9500 

1 1 .8 

4500 

7.9 

10000 

12.0 

5000 

6.4 





TABLE 3. Solar FJare Production 


Solar 

Flare 

cp_( >30 Mev) 

p -2 
p cm 

P 0 

Mv 

Yearly Average Production Rate 

r 1 4 -2 c -1 

C cm sec 

Normal 

Cutoff 

20% Normal 
Cutoff 

1956 

Feb 

23 

6.5x10® 

325 

0.86 

2.33 

1957 

Jan 

20 

3x10® 

60 

0.-306 

0.015 

1958 

Mar 

23 

2x10® 

55 

0.004 

0.009 


Jul 

7 

3x10® 

55 

0.005 

0.013 

1959 

May 

10 

7x10® 

60 

0.013 

0.034 


Jul 

10 

8.8x10® 

90 

0.039 

0.098 


Jul 

14 

1 . 1 x 1 0 9 

70 

0.030 

0.072 


Jul 

16 

8.1x10® 

110 

0.055 

0.14 

I960 

Nov 

12 

1 .4x1 0 9 

145. 

0.18 

0.45 


Nov 

15 

5.2x10® 

135 

0.056 

0.14 

1961 

Jul 

12 

1 .0x10® 

50 

0.002 

0.004 


Jul 

18 

2.1x10® 

135 

0.023 

0.057 







TABLE 4. 



Activity Increases from Supernova Gamma Ray Bursts 


(Konstantinov and Kocharov, 1967) 


Supernova Date A.D. Distance % Increase in 

( kpc) for W = 10 49 


Crab Nebula 

1054 

1 .1-1 .7 

0. 2-0.4 

Tycho Brahe 1 s 

1572 

0.36-3.3 

0.05-4 

Kepi er ‘ s 

1604 

1-9.9 

0.005-0.5 

Cassiopeia A 

1700 

3.4 

0.05 


Activity 
ergs . 
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Fi gure 1 . 
Figure 2. 


Fi gure 3 . 


Figure 4. 


Figure 5 . 


Figure 6. 


Radiocarbon yield as a function of incident particle 
rigidity for protons, alpha particles and CNO nuclei. 
Global average radiocarbon production as a function 
of the modulating parameters, n, or of the counting 
rates of the Mt. Washington neutron monitor and the 
Cheltenham ion chamber. Solar minimum, 1965, and 
solar maximum, 1967, correspond to n * 0 and 
r\ = - 2GV, respectively. 

Global average radiocarbon production as a function 
the effective geomagnetic dipole moment for various 
values of the modulating parameter, t). 

Variation of the radiocarbon excess in the upper 
reservoir as a function of time for the magnetic 
data given In Table 2. The measurements are normalized 
to the activity of 19th century wood. The parameter, 

R Q , i s the ratio between the present production and 
decay in the total reservoir. 

Global average radiocarbon production from solar flare 

particles as a function of P Q , defined in equation (19), 

o 

and normalized to 1 proton per cm sec with energy 
greater than 30 Mev. 

Variation of radiocarbon production over solar cycle 
19 due to solar flares and modulation by the inter- 
planetary field. The negative correlation of the 
radiocarbon product ion wi th solar activity, resulting 
from cosmic- ray modulation is obscured by enhanced 
particle fluxes from flares. 




Figure 7. 


Figure 8. 


Global average radiocarbon production from a super- 
nova source as a function time, normalized to the time 
to maximum at the earth ( t ) of relativistic particles 
from the supernova. 

Global average radiocarbon production as a function 
of time for cosmic rays from a supernova explosion 
superimposed on a constant cosmic- ray background flux. 
For a given distance, age, and diffusion mean free 
path the energy output of the supernova and the 
relative contribution of the background are such that 
during the last 8000 years the radiocarbon production 
decreased by 5% or 1% for the two curves, respectively. 
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